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We propose a new model based on cellular automata technique to simulate the behavior of Brow- 
nian 's particles with restrictions. In our model each particle will moves randomly for each time- step 
in each of his possible directions for a two-dimensional network. The movement of each particle 
is the result of collisions with other particles and this interaction determines the movement in one 
direction; they can only occupy one cell at time. We take only a representative percentage of fluid 
particles. We calculate the square root of the arithmetic mean of particles displacements (SRAMPD) 
in order to make a comparison with the diffusion coefficient in an adiabatic fluid. We also observe 
the behavior of the (SRAMPD) when we change the network size of the simulation for a fixed 
concentration and when we change the particles concentration for a fixed network size. 

PACS numbers: 87.55.kd, 05.40.Jc, 46.65.+g 



1. INTRODUCTION 

Brownian motion is the random movement for some 
nanoscopic particles in a fluid. It is named in honor 
of Robert Brown who described it in 1827. The ran- 
dom movement of these particles is the result of constant 
bombardment of his surface by fluid molecules under a 
thermal agitation. Atomic bombardment at this scale is 
not always completely uniform and have large statistical 
variations. The pressure on one side can vary, causing 
the movement observed [4]. 

The mathematical description of the phenomenon was 
developed by Albert Einstein [1]. Einstein found a way 
to confirm the atomic nature of matter observing the re- 
lationships between the macroscopic diffusion coefficients 
D and the atomic properties of matter. This relationship 
is: 



D = {RT)/{Na6r]7T) 



(1) 



Since R is the gas constant, Na Avogadro number, T is 
temperature in Kelvin, r] the viscosity, "a" the radius of 
the Brownian particles and D is the diffusion coefficient 
of the material suspended in the liquid [1,2]. The theory of 
Brownian motion was developed in order to describe the 
dynamic behavior of particles whose mass and size are 
much larger than the rest in the medium in which they 
are. Einstein was succeeded in proving that the average 
movement of the brownians particles in one direction is 
an expression such as: 



A = ^/l2Dt) 



(2) 
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In this work we present a new simple method to simu- 
late the diffusive behavior of particles that interact in a 
fluid, we control the statistical characteristics of the par- 
ticles movement by simples integers sums. We observed 
only the Brownians particles in the simulation. Each par- 
ticle is localized in a two-dimensional network and they 
have for each time-step one movement direction. This 
direction will be chosen randomly by our algorithm. 

The simulation represents a percentage of total parti- 
cles immersed in fluid, showing the phenomenon of self- 
diffusion. We perform calculations for the (SRAMPD) 
in order to compare with the experimental values ob- 
tained for the diffusion coefficient in [1]. We calculate 
the standard deviation (SD) of the values obtained for 
the (SRAMPD) in 50 simulations as a measure of the 
dispersion of the simulated data. Also we make simula- 
tions changing the concentration of particles fixing the 
network size, and changing the size of the network for a 
fixed concentration. In both cases we calculate the be- 
havior of the (SRAMPD) and the (SD). 



2. SIMULATION TECHNIQUE 

Our model use two networks represented as two matrix 
(Mi and Me). Each number inside on the matrix cell's in 
"Mi" represent the position in "Me" matrix. This num- 
ber may change for each time step and determine the par- 
ticle position in the spatial matrix Me. The supra-index 
of each cell in Mi represents the name of each particle in 
Me. (see Figure 1). This method will be called "Reticular 
Matrix Mapping" or MMR for his Spanish acronym. 

In each time step each particle will move randomly in 
one of the possible directions. Each particle will have 
4 possible movement directions for the two-dimensional 
case and 5 if we take the possibility of non- movement. To 
perform this movement, for each particle in Mi one may 
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FIG. 1: Example: Particle 1 (position 1x1) displayed in the 
cell as the number "4" this number represent the position in 
Me (1x4). The particle 3of Mi it is in the position 16 and so 
on for all other particles. The network size it's represented by 
N3xN4 = 5x5=25, and the numbers of particles by NlxN2 = 
3x3=9.. 







k 


± 


3 


4 


5 


i 


2 


3 


Me- 




7 






lO 


4 


5 




\\ 


12 


1^ 


14 


IS 








l> 




\^ 






8 / 












ST— 




S3 1 


J4- ■ 





13+5=18 



1 




6^ 


4 


5 


I a 


7 


8 


9 



1 


2 


3 


4 


5 


6 


7 


8 


9 


10 




12 


13 


14 


15 






\^ 


19 


2b 




22 


23 


24 


25 



FIG. 2: Example of particle 6 displacement from position 13 
to position 18. 



sum a integer, for each time-step. This sum may possible 
the particles movement to other position, to move any 
particle to the rigth side we need to add (+1) at the 
corresponding Mi cell, this will move the particle to the 
right side in Me matrix. If we sum (-1) we will move the 
particles to the left side, (-N3) for up, (+N4) for down 
and (+0) for non-movement. Nl and N2 are the matrix 
components for Mi and N3, N4 the matrix components 
in Me. 

As an example: We take the particle 6 in position 13, 
if we add (+5) to the particle 6 we will move the particle 
from position 13 to 18 in "Me" (see Figure 2). 

Our particle systems have restrictions to moving up, 
down, right, left and no movement. We can assign a per- 
centage of movement for each direction. Every time- step 
each particle will have the chance to move in a direc- 
tion that will be chosen randomly. We assign a weight 
of probability for each of his possible directions. In this 
case for our propose, all direction must have the same 
probability weight to correspond a random movement in 
all directions like in a Brownian movement or a 2 dimen- 
sional random walk. (20% up, down, right, left and non 
moving) . 

This method allows us to stadistastically control the 



FIG. 3: Steps 1, 2, 3, 4 for a simple system simulation 
with Me=5x5=25cells, Mi= 1x1=1 particle with 5 directions 
of movement (up and down left and right zero) with 20% 
probability of each movement direction with closed borders. 




FIG. 4: Simulation in 1000 time steps. Me have 
150xl50=22500cells, and Mi=50x50=2500 particles. We use 
the probability weights 30% down, 20% up, 20% right, 20% 
left and 10% zero movements. For many steps in time, almost 
all the particles are located at the bottom of the simulation. 
The first figure from left to rigth correspond at the first time- 
step, the second the lOOstep, 500 and lOOOstep. 



particles movement. In the Figure 4 we show snap-shots 
for particles which moves randomly with different weight 
porcentage. In this case, the particles will move prefer- 
entially to the bottom part, like high density particles in 
a classic fluid, subject only to gravity forces, i.e. A more 
higher probability to move in a "down" direction. The 
simulation showed in the figure 4 correspond with weight 
porcentage probabilities: down 30%, up 20%, left 20%, 
right 20% and non-movement of 10% in 1000 steps sim- 
ulation. We can appreciate a clear trend for all particles 
to stay in the bottom part of the picture. 

A mechanism in the code prevents particles with the 
same number in the Mi matrix, avoiding particles in the 
same position. For the boundary conditions are estab- 
lished rules on the border cells to determine the type 
of boundaries that we want. In our case, we choice as 
boundary conditions rigid-walls enclosing all particles. 



3. SQUARE ROOT OF THE ARITHMETIC 
MEAN OF PARTICLES DISPLACEMENTS 



In order to compute (SRAMPD) we take the initial and 
final position for all particles in the network. The number 
of particles is given by " n" sub-index. The initial position 
are represented by " i" sub-index and the final position by 
the "f sub-index . Where n = NlxN2 is the number of 
particles in Mi. Using the familiar expression for the two- 
dimensional distance between two points; summing and 
dividing over all particles we obtain the (SRAMPD) in 
the simulation: 
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FIG. 5: Values of 7 obtained in 50 simulations with the same parameters Mi 30x30 = 900 particles, Me 100x100=10.000 cells 
in 200 time-steps. The probability for each movement direction was 20% for each one: Up, down, right, left and not moving. 
The mean of 7 in the computational experiments carried was 7=10.3159 cell units and the standard deviation was 0.0152 in 
the same units. 
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FIG. 6: 7 Values for different network sizes with the same particles concentration. The particles concentration was fixed at 
20.25%. If Me = 6400 = 80x80particles, and Mi=36x36=1296 particles the concentration will be 20.25%. Each data obtained 
for 7 was the average of 50 simulations. Each simulation was performed in 200 time-steps with the same characteristics. The 
probability for each motion direction was 20% up, down, left and right and non movement. 



(3) 



To comparate (SRAMPD) in the simulation 7, whit 
the experimental (SRAMPD) A we need to establish a 
scale parameter a associating in this, the simulation scale 
with the real scale. This expression will be given by the 
product between the (SRAMPD) in the simulation 7 and 
the scale factor a: 



A = (77 



If N is equal to 6x10^ in accordance with the kinetic 
theory of gases, water at IT^'C, k = 1.35x10"^^ J/iC, the 
particles diameter 0,001mm we have that A in the x di- 
rection is 0.8/i=8xl0~^cm and the displacement of a half 
minutes will be 0.6/i [1]. This value is independent of 
the number of experiments made it and its a invariable 
value. 



If we take the (SRAMPD) simulation value: 
7=10.3159 and divide between the experimental value 
of A = 6/i in [1]. We can calculate the value of a using 



(4) (4) a = 0,5816xlO-V. 



Figures 6 and 7 show the behavior of 7 when we vary 



the network size maintaining the same particles concen- 
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FIG. 8: Graphic for 7 values varying the concentration of particles and fixing the size of the network Me= 100x100= 10. OOOcells. 
Each measure is the average obtained from 50 simulations for 7, each simulation was performed under the same conditions 20% 
probability for each movement direction: up, down, left and right does not move in 200 time-steps. 
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FIG. 7: Standard deviation (SD) for the data obtained in the figure 6 for all of the 50 simulations in 7. These values represent 
a measure of data dispersion in the simulation. 



tration 20.25% for each simulation. This concentration 
is a measure between the total number of cells in Me 
and the number of particles in Mi. Each data obtained 
come from a average of 50 simulations performed with 
the same characteristics for 7; 200 time-steps and the 
movement probability for each direction was 20%: up, 
down left and right and non movement. The results in 
Figure 6 shows a similar behavior to Hardy and Pomeau 
[7] they measured the average time of free flight in a 
network, fixing the concentration and increasing size of 



the network. His results are comparable with our re- 
sult however our method propose a direct way to change 
the movement statistical characteristic in the simulation 
and a new simulation technique. Figure 7 shows how the 
scattering of data decreased as the network gets bigger, 
larger networks produce less dispersion in data. The dis- 
persion value in the data are small if we compare whit 
the size of the network and the values of (SRAMPD). 
The minimum was approximately 0.02 cell units and the 
maximum 0.14 cells units. 
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Graphs 8 show the behavior of the (SRAMPD) when 
we varying the concentration of particles in the simula- 
tions and fix the network size. We obtained data for 18 
different concentrations ranged from 1% to 90% of oc- 
cupancy. Each point corresponded to the average value 
of 7 in 50 simulations; each simulation was performed 



in 200 time-steps with the same movement probability: 
up, down, left right and non movement 20%. The net- 
work size for all cases was 100x100 cells or 10.000 possible 
positions. We see as the value of (SRAMPD) decreases 
when we increasing the particle concentration. That is 
attributed to the fact that the movement of particles 
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FIG. 9: (SD) for the data obtained in the figure 8 each point correspond to the dispersion of 50 simulations in 7. These values 
represent a measure of the dispersion data in the experiments. 



depends of the numbers of available possibilities in his 
neighbors, if we have a more higher numbers of particles, 
the possibilities will be decrease changing the particles 
mobility. The values of (SD) remained in a relatively 
small variations if we compared with the size of the net- 
work and the simulation magnitud of (SRAMPD). 

4. CONCLUSIONS 

The properties on invariability obtained by Einstein in 
[1] for the compute of the (SRAMPD) and the respective 
comparison with the diffusion coefficient in a fluid, give 
us a prove of the randomness characteristic of the Brow- 
nian movement phenomenon. This is the same principle 
used in random walks in 2D. No matter how many times 
does the experiment we always get the same results. We 
represent this situation with our new method MMR and 
we obtained the same results for many experiment. The 
graphics 6 show us this result with a little dispersion 



proved by the smalls values of (SD). For 50 experiments 
we obtain a dispersion of 0.0152 cells units and the mean 
in 50 simulations was 7=10.3159cells units. In order to 
compare with real experiments showed in [1] we propose 
a scale parameter a. We perform simulations changing 
the size of the network and fixing the concentration, and 
fixing the network size and changing the concentration. 
In the first case we observe an increase of 7 and a soft 
decline in the dispersion. In the second case we observe 
a decline of 7 caused by the increasing number of par- 
ticles. For too many particles they lose possibilities to 
occupy positions and mobility in general decrease. The 
method allows us to implement a so higher complex sit- 
uation, and this work represents only ones of the more 
simples experiment. We will go to per forme simulations 
in: the particles movement having preferences to move in 
one special direction, simulations taking in a count the 
collisions between particles, 3D networks, fluid particles 
moving in a porous media, and so other situations. 
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